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Abstract 

5h The fusion coefficients from 50(3) to 5*0(4) play a key role in the definition of spin foam models 

for the dynamics in Loop Quantum Gravity. In this paper we give a simple analytic formula of 
the EPRL fusion coefficients. We study the large spin asymptotics and show that they map SO{3) 
, semiclassical intertwiners into SU{2)l x SU{2)r semiclassical intertwiners. This non-trivial property 

^ ' opens the possibility for an analysis of the semiclassical behavior of the model. 

00 " 

■ 1 Introduction 

cn 

. The recent construction of a class of spinfoam models [U |2l El SI |5] compatible with loop quantum 

C3 ' gravity (LQG) [6, 7, 8, 9J has opened the possibility of consistently defining the LQG dynamics using 

00 , spinfoam techniques [lOj |TTJ [T^ [T3] • In this paper we focus on the Engle-Pereira-Rovelli-Livine (EPRL) 

spinfoam model for Riemannian gravity introduced in [3]. For given Immirzi parameter 7, the vertex 
amplitude is defined as follows: it is a function of five 5*0(3) intertwiners ia and ten spins jat (with 
I a, 6 = 1, . . , 5 and a < b) given by 

The functions {15j}jv are normalized 15j-symbols, namely the contraction of five normalized 4-valent 
SU{2) intertwiners according to the pattern of a 4-simplex, and the fll^a are fusion coefficients from 
50(3) to SU{2)l X SU{2)ji introduced in [3] and defined below. Such coefficients play a key role in 
the definition of the model. Indeed the model differs from the one introduced by Barrett and Crane 
[13] only for the structure of these coefficients. In this paper we study the large spin asymptotics of the 
EPRL fusion coefficients. 

A careful analysis of the asymptotics of fusion coefficients is a step needed for the study of the 
semiclassical properties of the model. In fact, we have already used the results that we present in this 
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paper in order to understand the features of the wavepacket evolution. The propagation of semiclassical 
wavepackets was introduced in [15j as a new way to test the semiclassical limit of a spinfoam model. 
A spinfoam model has a good semiclassical behavior if semiclassical wavepackets (peaked on a classical 

3- geometry) follow the trajectories predicted by the classical equations of motion. In [15] this new 
technique was implemented in the EPR flipped vertex model to study the propagation of intertwiner 
wavepackets. In |16j we developed a more efficient numerical algorithm (using techniques similar to 
|17j ) and applied the asymptotic analysis presented here. This kind of study is, in the general context 
of spinfoam models, complementary to the semiclassical analysis based on the calculation of n-point 
functions. 

In [18| [19] , a strategy for recovering graviton correlations from a background-independent theory was 
introduced. The idea was tested on the Barrett-Crane model at the "single-vertex" level. At this level, 
correlations of geometric operators can be checked against perturbative Regge-calculus with a single 

4- simplex [3D]. Given the fact that the Barrett-Crane model gives trivial dynamics to intertwiners, 
the analysis was restricted to the spin degrees of freedom - namely to area correlations only. On the 
other hand, the new models are consistent with the LQG kinematics and allow the computations of 
semiclassical correlations of geometric observables as the area, the angle, the volume or the length 
[2T] [22l [23l [24l [25l [26] . At the single- vertex level, the semiclassical correlations for two local geometric 
operators Oi, O2 are simply given by 

where W^(ja6,*a) is the vertex-amplitude introduced in ([1]) and ^'g(ja6,*a) is a boundary semiclassical 
state peaked on a configuration q of the intrinsic and the extrinsic geometry of the boundary of a region 
of space-time. The appropriate dependence on spins and intertwiners of the state 5'g(jab, ia) is discussed 
in [27l [28] and uses the semiclassical tetrahedron state of [29] . Moreover, in order to guarantee that 
the appropriate correlations are present, in [571 dH] a specific form of the large spin asymptotics for 
the vertex amplitude was conjectured (see [3D]). In order to show that the EPRL vertex amplitude 
satisfies this conjecture, an analysis of the asymptotics of the fusion coefficients is needed. The region 
of parameter space of interest is large spins jab and intertwiners ia of the same order of magnitude of 
the spins. As a result, the fusion coefficients for the node a, f^l^nijab), can be seen as a function of the 

two bare variables , i§, of the fluctuation of the intertwiner ia and of the fluctuation of the four spins 
jab- In this paper we focus on this analysis. For different approaches to the semiclassical limit, see [3T] 
and [31]. 

The paper is organized as follows: in section [D we show a simple analytic expression for the EPRL 
fusion coefficients; in section[3]we use this expression for the analysis of the asymptotics of the coefficients 
in the region of parameter space of interest; in section[4|we show that the fusion coefficients map 5*0(3) 
semiclassical intertwiners into SU{2)l x SU{2)fi semiclassical intertwiners. We conclude discussing the 
relevance of this result for the analysis of the semiclassical behavior of the model. In the appendix we 
collect some useful formula involving Wigner coefficients. 



2 Analytical expression for the fusion coefficients 

The fusion coefficients provide a map from four-valent iS'0(3) intertwiners to four-valent SO{A) inter- 
twiners. They can be defined in terms of contractions of SU{2) 3j-symbols. In the following we use a 
planar diagrammatic notation for SU{2) recoupling theory [33]. We represent the SU{2) Wigner metric 
and the SU (2) three- valent intertwiner respectively by an oriented line and by a node with three links 
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oriented counter-clockwise^. A four-valent iS'0(3) intertwiner \i) can be represented in terms of the 
recoupling basis as 



,72 .j3\ 

where a dashed hne has been used to denote the virtual link associated to the coupling channel. Similarly 
a four-valent 5*0(4) intertwiner can be represented in terms of an SU{2)l x SU{2)ii basis as |«L)|i_R). 

Using this diagrammatic notation, the EPRL fusion coefficients for given Immirzi parameter 7 are 
given by 

fU^{3i.j2,3z,jA) =(-l)^'^-^'^+^'^-^y(2z + l){2iL + l){2iR + l)K^,{2u + 1) x (4) 



|l-7|j 




1 — rlia (i+7)j4 
I 2 I 2 I 




These coefficients define a map 
/ : Inv[Hj^ (g) . . . ( 



H 



Inv\H, {i+2)2x\ 
y 2 ' 2 ) 



H, ii-t|j4 a+2W4 J 
\ 2 ^ 2 I 



(5) 



from 50(3) to 5*0(4) intertwiners. Using the identity 



I I 



(6) 



where the shaded rectangles represent arbitrary closed graphs, we have that the diagram in ([4]) can be 
written as the product of two terms 

/^.«(J1,J2,J3,J4) = (-l)^-i-^-^+J3-.4^(2, + 1)(2,^ + l)(2z^ + l)n„(2j„ + 1) (7,\,Jjl,j2) gL..03,j4) 

(7) 

where ql^i^ is given by the following 9j-synibol 



*„(ii>j2) 




— 1 — r|j2 



^Ji ^L ^^2 



(8) 



Ji 



J2 



ih ^ iR 



minus sign in place of the + will be used to indicate clockwise orientation of the links. 
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From the form of qlj^i^^ we can read a number of properties of the fusion coefficients. First of aU, 
the diagram in expression ([S]) displays a node with three links labelled ijiLjiR- This corresponds to 
a triangular inequality between the intertwiners which is not evident from formula As a 

result we have that the fusion coefficients vanish outside the domain 

|«L - ifll <i <iL + iR ■ (9) 

Moreover in the monochromatic case, ji = j2 = J3 = ji, we have that the fusion coefficients are non- 
negative (as follows from ([7])) and, for il + + * odd, they vanish (because the first and the third 
column in the 9j-symbol are identical). 

As discussed in [3J[S], the fact that the spins labeling the links in ^ have to be half-integers imposes a 
quantization condition on the Immirzi parameter 7. In particular 7 has to be rational and a restriction 
on spins may be present. Such restrictions are absent in the Lorentzian case. Now notice that for 
< 7 < 1 we have that -I- ^"'""^^ — 1 , while for 7 > 1 we have that — ■'^^-p^ — 1 (with the limiting 
case 7 = 1 corresponding to a selfdual connection). As a result, in the first and the third column of 
the 9j-symbol in the third entry is either the sum or the difference of the first two. In both cases 
the 9j-symbol admits a simple expression in terms of a product of factorials and of a 3j-symbol (see 
appendix Using this result we have that, for < 7 < 1, the coefficient ql^i^^iii, j2) can be written 
as 



9L.Jji,j2) = (-i)— (^^^^^ -(,,-,2)j ^--(^^'^^^ 

^it^ (ii'.?2) given by 



A,z«Ui:J2)- \l ^o.. +1)1(2^2 + 1)! ^^^> 




(|l-7b-i)!(|l-7b2)! 



+ _ i^) I _,_ Mi + + 1) ! 



((l+7)ji)!((l + 7)j2)! 



A similar result is available for 7 > 1. The Wigner 3j'-symbol in expression (|10p displays explicitly 
the triangle inequality ([9]) among the intertwiners. Notice that the expression simplifies further in the 
monochromatic case as we have a 3j-symbol with vanishing magnetic indices. 

The fact that the fusion coefficients ^ admit an analytic expression which is so simple is certainly 
remarkable. The algebraic expression ((7l). p0|) . (|lip involves no sum over magnetic indices. On the other 
hand, expression (jj]) involves ten 3j-symbols (one for each node in the graph) and naively fifteen sums 
over magnetic indices (one for each link). In the following we will use this expression as starting point 
for our asymptotic analysis. 



3 Asymptotic analysis 

The new analytic formula (Tf]! . (fTO|) . (fTT|) is well suited for studying the behavior of the EPRL fusion 
coefficients in different asymptotic regions of parameter space. In this paper we focus on the region 
of interest in the analysis of semiclassical correlations as discussed in the introduction. This region 
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is identified as follows: let us introduce a large spin jo and a large intertwiner (i.e. virtual spin in a 
coupling channel) iq; let us also fix the ratio between io and jo to be of order one - in particular we will 
take iQ = -^joi then we assume that 

• the spins ji, j2, js, ji, are restricted to be of the form je = jo + ^je with the fluctuation Sje small 
with respect to the background value Jq. More precisely we require that the relative fluctuation 
^ is of order o{l/y/%); 

• the SO{3) intertwiner i is restricted to be of the form i = io + Si with the relative fluctuation ^ 
of order o(l/Vjo); 

• the intertwiners for SU{2)l and SU{2)fi are studied in the region close to the background values 
^0 _ jp and i^j = ^^-^io- We study the dependence of the fusion coefficients on the fluctuations 
of these background values assuming that the relative fluctuations Sil/iq and Sir/iq are of order 
o(l/\^). 

A detailed motivation for these assumptions is provided in section |4l Here we notice that, both for 
< 7 < 1 and for 7 > 1, the background value of the intertwiners «l, za, i, saturate one of the 
two triangular inequalities ([9|). As a result, we have that the fusion coefficients vanish unless the 
perturbations on the background satisfy the following inequality 

5i<5iL+5iR 0<7<1 (12) 
5iR<5i + 5iL 7>1. (13) 

In order to derive the asymptotics of the EPRL fusion coefficients in this region of parameter space we 
need to analyze both the asymptotics of the 3j-symbol in PH)) and of the coefficients j2). This 

is done in the following two subsections. 



3.1 Asymptotics of 3j-symbols 

The behavior of the 3j-symbol appearing in equation ((T0|) in the asymptotic region described above is 
given by Ponzano-Regge asymptotic expression (equation 2.6 in [34]; see also appendix [B|) : 



^L 



|l-7l(3l-j2) (l+7)(jl-32) 
2 2 



-(ji ~h) 



(14) 



(-1) 



V27rA 



■ cos 



(< 



|l-7l(jl-32) ^ 



(l+7)(jl-j2) ^ 



4/ 
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The quantities A, 9l, Or, 0, 0-, admit a simple geometrical representation: let us consider a triangle 
with sides of length i_L + |, i_R + ^, i + \ embedded in 3d Euclidean space as shown below 




In the figure the height of the three vertices of the triangle with respect to a plane are given; this fixes the 
orientation of the triangle and forms an orthogonal prism with triangular base. The quantity A is the 
area of the base of the prism (shaded in picture). The quantities Ol, 6*^, 6 are dihedral angles between 
the faces of the prism which intersect at the sides il, ir, i of the triangle. The quantities (j)^ are 
dihedral angles between the faces of the prism which share the side of length h + \1 — — j2)/2 and 
the side of length h -\- (1 + 7) (ji — j2)/2, respectively. For explicit expressions we refer to the appendix. 

In the monochromatic case, ji = j2 j we have that the triangle is parallel to the plane and the formula 
simplifies a lot; in particular we have that the area A of the base of the prism is simply given by Heron 
formula in terms of il, iR, i only, and the dihedral angles 9l, Or, are all equal to 7r/2. As a result the 
asymptotics is given by 

1 l + (16) 




2 

Notice that the sum il + iR + i is required to be integer and that the asymptotic expression vanishes 
if the sum is odd and is real if the sum is even. Now, the background configuration ol il, ir and i we 
are interested in corresponds to a triangle which is close to be degenerate to a segment. This is due to 
the fact that ^^^-j^io + ^"^^^^ io = io for < 7 < 1, and ^^^yi^io — ^'^^ ^^ io = »o for 7 > 1- In fact the 
triangle is not degenerate as an offset ^ is present in the length of its edges. As a result the area of 

3/2 

this almost-degenerate triangle is non-zero and scales as ig for large iq. When we take into account 
allowed perturbations of the edge-lengths of the triangle we find 

[\VT^il'^.JTV2{5nT5iI^)+o(i^^'^)) 0<7<1 

A={ (17) 

[jy/l^^T (Vl + 2('^* + 5iL - Sir) + o(*o '/")) 7>1 ■ 

This formula holds both when the respective sums SiL + ^iR — Si and Si + SiL — Sin vanish and when they 
are positive and at most of order 0(\/zo). As a result we have that, when SiL + Sin — Si, or Si + Si^^Sin 
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respectively, is even the perturbative asymptotics of the square of the Sj-symbol is 

2 

I - (18) 



2 "U ' "-I. 2 





e{SiL+SiR-di) 0<7<1 



i /1 ^wlC , ■ , ^(^» + - 7>1 ■ 

(19) 

The theta functions implement the triangular inequality on the fluctuations. In the more general case 
when ji — j2 is non-zero but small with respect to the size of the triangle, we have that the fluctuation in 
dje can be treated perturbatively and, to leading order, the asymptotic expression remains unchanged. 

3.2 Gaussians from factorials 

In this subsection we study the asymptotics of the function A\^^^{ji, j2) which, for < 7 < 1, is given 
by expression . The proof in the case 7 > 1 goes the same way. In the asymptotic region of interest 
all the factorials in (jlip have large argument, therefore Stirling's asymptotic expansion can be used: 

N 

j„! = e+^o(l°g^o - 1) {l + Y^ a„jV" + 0(jo"^^+^^)) for aU iV > 0, (20) 

where a„ are coefficients which can be computed; for instance ai ^ j^. The formula we need is a 
perturbative expansion of the factorial of (1 + C)jo when the parameter ^ is of order o{\/^/Jq). We have 
that 

CXD 

((1 + Ok) ! = exp ( + io(log JO - 1) + Oo log JO + Jo CkS!') x (21) 

fc=i 

N M 

X (1 + E E ""^"^o"" ^" + OU"^"^^^^'^)) (22) 



where the coefficients hm and Ck can be computed explicitly. We find that the function ^i^i^(ji, J2) has 
the following asymptotic behavior 

^Sn,,,- c™^,,- (jo + <5ji,jo + <5j2) ^ Ao(jo)e-^(^*-^*-^^^^^-^^-^) (23) 
where Ao(jo) is the function evaluated at the background values and H{5iL, Siji, 5i, Sji, 5j2) is given by 
H{6iL,SiR, Si, Sji,Sj2) = ^(arcsinhVS) {6iL + Sir - 5i)+ (24) 
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2 |l-7|io 2 (1+7)10 4 io 

l5iL+6iR-6i 1 

-(oji + '5j2) + 0{—j=) 



2 io Vjo 
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for < 7 < 1, while for 7 > 1 it is given by 



H{SiL, SiR, Si, Sji, Sj2) = - (arcsinhA/S) {Si + SiL — SiB)+ (25) 



2 |l-7|io 2 (l+7)«o 4 

1 (5i + (5iL - (5ifl 1 

-[Sji +Sj2) + 0{—=) 



2 io Vjo' 



3.3 Perturbative asymptotics of the fusion coefficients 

Collecting the results of the previous two subsections we find for the fusion coefficients the asymptotic 
formula 

(i+.).o Oo + 6j,) ^ /oOo) -^===L==== d{SiL + Sir - Sz) x (26) 

2 +'^'^^ 2 yJl + 2{SlL+SlR- 5i) 

X exp ( — arcsinh(\/3) (Jil + — x 



X exp ( - v3-- V3-; — + — : — ) X 

^ II-7U0 (l + 7)«o 2 iQ ' 



,lSiL + Sir- Si , X - \\ 

X exp (- : {Sji + Sj2 + Sj3 + Sj4)) 

for < 7 < 1, and 

fi"-'k (i+.).c Oo + 6je) ^ foUo) -j===l==== e{Si + SlL - Sir) X (27) 

X exp ( — arcsinh(\/3) (Si + — SiR)) x 

(<5*i)2 /- (^z;,)2 V3(<5*)% 



X exp 



l-7Un 



II-7N0 (l+7)«o 2 io 



X exp (- : (dji + d:/2 + dj3 + <3j4)j 

for 7 > 1, where /o(jo) is the value of the fusion coefficients at the background configuration. As we will 
show in next section, this asymptotic expression has an appealing geometrical interpretation and plays 
a key role in the connection between the semiclassical behavior of the spin foam vertex and simplicial 
geometries. 



4 Semiclassical behavior 

In [15] the propagation of boundary wave packets was introduced as a way to test the semiclassical 
behavior of a spinfoam model. In particular, the authors considered an "initial" state made by the 
product of four intcrtwiner wavepackets; this state has the geometrical interpretation of four semiclassical 
regular tetrahedra in the boundary of a 4-simplex of linear size of order Vjo- Then this state was evolved 
(numerically) by contraction with the flipped vertex amplitude to give the "flnal" state, which in turn is 
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an intertwiner wavepacket. While in '15] only very small jo's were considered, in [16j we make the same 
calculation for higher spins both numerically and semi-analitically, and the results are clear: the "final" 
state is a semiclassical regular tetrahedron with the same size as the incoming ones. This is exactly 
what we expect from the classical equations of motion. 
The evolution is defined by 

X! Wijo,h,---,i5)ip{h,jo)---ip{i5,jo) = H'i5,jo), (28) 

ii...i6 

where 

i'ihJo) = C{jo) exp ( - ^ll^l^ + _ (29) 

is a semiclassical 50(3) intertwiner (actually its components in the base |i)), or a semiclassical tetrahe- 
dron, in the equilateral configuration, with C{jo) a normalization constant, and W{jo,ii, ... ,15) is the 
vertex (jT]) with 7 = evaluated in the homogeneous spin configuration (the ten spins equal to jo). In 
(|28p . if we want to make the sum over intertwiners, for fixed jo, then we have to evaluate the function 
g defined as follows: 

giiL,iR,jo) = f'LiRijo)i^{hjo) ■ (30) 

i 

The values of g are the components of an 5*0(4) intertwiner in the basis |jL)|j_R), where is an inter- 
twiner between four SU{2) irreducible representations of spin Jq = -^^-j^joj SLud is an intertwiner 
between representations of spin = ^^jo- 

We show that EPRL fusion coefficients map 50(3) semiclassical intertwiners into SU{2)l x SU{2)ji 
semiclassical intertwiners. The sum over the intertwiner i of the fusion coefficients times the semiclas- 
sical state can be computed explicitly at leading order in a stationary phase approximation, using the 
asymptotic formula (pS))(|27p . The result is 

I1-7I 



E M ^i^^Jo) « ao foUo) Cijo) X exp ( - ^3^^^^-^^-^ ± i|(*i - M x (31) 



(1+7) 
2 

(1+7) 

where 



X exp ( - V3 + .-{tn - ^^o) 



p — arcsinh( V^) 

ao = E .i,2fc " ' 

the plus-minus signs both in ([3T|) and ((32)) refer to the two cases 7 < 1 (upper sign) and 7 > 1 (lower 
sign). The r.h.s. of (j3ip . besides being a very simple formula for the asymptotical action of the map / 
on a semiclassical intertwiner, is asymptotically invariant under change of pairing of the virtual spins ii, 
and ijf (up to a normalization iV). Recalling that the change of pairing is made by means of 6j-symbols, 
we have 




EE v/dim(.L)dim(*H)(-l)-+'=^+-+^- \ lAi |x (33) 

g{iL,iR) ^ N(jo) g(kL,kR,jQ) . 
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(a) (b) 

Figure 1: (a) Interpolated plot of the modulus of g{iL, iR,ja) for = 20 and 7 = computed using the 
exact formula of the fusion coefficients, (b) Top view of the imaginary part. 



This result holds because each of the two exponentials in (|3T|) is of the form 

exp(-f^^±i|(fc-M), (34) 

which is a semiclassical equilateral tetrahedron with area quantum numbers /cq; it follows that g is 
(asymptotically) an SO{A) semiclassical intertwiner. The formula (|31|1 can be checked against plots of 
the exact formula for large jo's; a particular case is provided in fi-gH] 

In addition, we can ask whether the inverse map has the same semiclassical property. Re- 
markably, the answer is positive: maps semiclassical 50(4) intertwiners into semiclassical SO{?>) 
intertwiners. The calculation, not reported here, involves error functions (because of the presence of the 
theta function) which have to be expanded to leading order in 1/jo- 

A final remark on our choice for the asymptotic region is needed. The goal we have in mind is 
to apply the asymptotic formula for the fusion coefhcients to the calculation of observables like ^ in 
the semiclassical regime. If the classical geometry q over which the boundary state is peaked is the 
geometry of the boundary of a regular 4-simplex, then the sums in ([2]) are dominated by spins of the 
form jab = jo + 5 jab and intertwiners of the form ia = io + Sia, with io — 2 jo / \/3, where the fluctuations 
must be such that the relative fluctuations Sj / jo, Si/ jo go to zero in the limit jo 00. More precisely, 
the fluctuations are usually chosen to be at most of order 0(\/jo)- This is exactly the region we study 
in this paper. As to the region in the {iLiiR) parameter space, the choice of the background values 
J^i-^ioi ^~Y^^o and the order of their fluctuations is made a posteriori both by numerical investigation 
and by the form of the asymptotic expansion. It is evident that the previous considerations hold in 
particular for the function g analyzed in this section. 

5 The case 7=1 

When 7 = 1 we have that j^ = ^"'"^^^ j = and we can read from the graph @ that the fusion coefhcients 
vanish unless = 0. Furthermore it is easy to see that for 7 = 1 the fusion coefficients vanish also 
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when iR is different from i. This can be seen, for instance, applying the identity 



V7 

dun I Z^"*^ 



to the graph Q with ~ Q. As a resuh, we have simply 

/'i,i„(il,j2,j3, j4) = k^fiSia.t (36) 

and the asymptotic analysis is trivial. The previous equation can be also considered as a normalization 
check; in fact, with the definition (jl]) for the fusion coefficients, the EPRL vertex amplitude ([1]) reduces 
for 7 = 1 to the usual 50(3) BF vertex amplitude. 

6 Conclusions 

We summarize our results and give some outlook in a few points. 

• We have shown a simple analytic formula for the LQG fusion coefficients, as defined in the EPRL 
spinfoam model. 

• We have given a large spin asymptotic formula for the coefficients; specifically, we made a pertur- 
bative asymptotic expansion around a background configuration dictated by the kind of boundary 
state considered. 

• The picture coming out from our analysis is promising: the fusion coefficients not only give non- 
trivial dynamics to intertwiners at the quantum level, but they seem to behave very well at 
semiclassical level, in fact they map semiclassical 50(3) tetrahedra into semiclassical 50(4) tetra- 
hedra. This is to us a highly non-trivial property which, in turn, makes the semiclassical analysis 
of dynamics less obscure. A first application of the asymptotic formula can be found in [16j . 

• Our analysis is a step needed for the study of the full asymptotic expansion of the EPRL vertex, 
which is part of our work in progress. 
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A Properties of 9j-symbols 



The 9j-symbol with two columns with third entry given by the sum of the first two can be written as 
a f c 

9 d \ = (^-i)f-9+a+b-{c+d) , - ) X (37) 



a + b h c + d 



/ (2a)!(2b)!(2c)!(2d)!(o + b + c + d- h)l{a + b + c + d + h + 1)\ 

Y (2a + 2b+ l)!(2c +2d+ l)!(a + c - f)\{a + c+f + + d- g)\{b + d + g + l)\ ' 

An analogous formula for the 9j-symbol with two columns with third entry given by the difference of 
the first two can be obtained from the formula above noting that 




(38) 



so we are in the previous case. 

The 3j-symbol with vanishing magnetic numbers has the simple expression 



i ' { ^-"-''^^ y.^ia + b-cy. V (£+|^)!(£^)!( °+''+^+i )! 



(39) 



These formula can be derived from [331 ISS] . 

B Regge asymptotic formula for 3j-symbols 

The asymptotic formula of 3j-symbols for large spins a, b, c and admitted magnetic numbers, i.e. ma 
nib + TTic — 0, given by G. Ponzano and T. Regge in [34J is 



a b c 
ma mb mc 



{-\Y+b-c+l /111 7r\ 

COS [{a + -)ea + {b+ -)eb + (c + -)ec + ma(f>a - mb(t)b + - j 



with 



arccos {2{a + \fmc + TOa((c + \f + (a + \f - (6 + \f)^ 
y ((a + \Y - ml) (4(c + \Y{a + \Y - ((c + i)^ + (a + \Y - (b + i)^)^) 



(41) 



arccos -- 
2 



1 (a + i)2 - (6 + 1)2 - (c + 1)2 - 2mfcmc 



^((fe+i)2-m2)((c+i)2-m2) 



A: 



\ 



det 

16 



V 





(aH 






- W ml 




(aH 









-\f-ml 


1 




-i)2-m2 (cH 


- W ml 







1 




1 


1 




1 


0/ 



(42) 



(43) 



and ^b, 0c, 06 are obtained by cyclic permutations of (a, 6, c). 
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